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Abstract — Singular Value Decomposition (S VD) is a powerful 
tool for multivariate analysis. However, independent computation 
of the SVD for each sample taken from a bandlimited matrix 
random process will result in singular value sample paths whose 
tangled evolution is not consistent with the structure of the under- 
lying random process. The solution to this problem is developed 
as follows: (i) a SVD with relaxed identification conditions is 
proposed, (ii) an approach is formulated for computing the SVD's 
of two adjacent matrices in the sample path with the objective 
of maximizing the correlation between corresponding singular 
vectors of the two matrices, and (iii) an efficient algorithm is given 
for untangling the singular value sample paths. The algorithm 
gives a unique solution conditioned on the seed matrix's SVD. Its 
effectiveness is demonstrated on bandlimited Gaussian random- 
matrix sample paths. Results are shown to be consistent with 
those predicted by random-matrix theory. A primary application 
of the algorithm is in multiple-antenna radio systems. The benefit 
promised by using SVD untangling in these systems is that the 
fading rate of the channel's SVD factors is greatly reduced so 
that the performance of channel estimation, channel feedback 
and channel prediction can be increased. 

I. Introduction 

Many observable processes in the natural world give rise 
to matrix random process. Examples can be found in multiple 
antenna radio communications [1] [2], gene expression [3] [4], 
and kinematics [5] among others. If these random process are 
bandlimited, then they may be sampled at uniform intervals 
and represented as discrete sample paths without loss of 
information. The Singular Value Decomposition (SVD) is a 
powerful tool for analyzing the multidimensional observations 
from these matrix random processes. However, computation 
of the SVD for each matrix-valued sample of a bandlimited 
random process results in singular value sample paths whose 
evolution is not consistent with the structure of the underlying 
random process. This is because the SVD computation is 
applied independently from sample to sample, during which 
strict phase and ordering identification conditions are imposed 
on the singular values. 

This work addresses the problem of how to compute the 
SVD's of a random matrix sample path in a way that preserves 
the covariance of the underlying random process. The solution 
to this problem is developed as follows: (i) a SVD with relaxed 
identification conditions is proposed, (ii) an approach is formu- 
lated for computing the SVD's of two adjacent matrices in the 



sample path with the objective of maximizing the correlation 
between corresponding singular vectors of the two matrices, 
and (iii) an efficient algorithm is given for untangling the 
singular value sample paths. 

Random matrix processes of the type arising in Multiple 
Input Multiple Output (MIMO) radio channels will be studied 
as an example throughout this paper. This is a useful case 
to study because it lends an intuitive physical interpretation 
of the SVD and because it is a relevant research area. Spatial 
multiplexing of information across antennas in an array allows 
MIMO transceivers to communicate information over indepen- 
dent parallel channels. The SVD of a MIMO channel matrix is 
used to provide the transmitter and receiver with the necessary 
beamforming coefficients for spatial multiplexing [6]. 

The remainder of the work is organized as follows. In 
Section II, the traditional definition of the SVD is revisited 
and MIMO beamforming is introduced. Section III details the 
approach for computing the SVD's of two adjacent matrices 
and the untangling algorithm is presented. Section IV demon- 
strates the effectiveness of the algorithm by presenting results 
from the untangling of a MIMO Gaussian random process. A 
summary of the major contributions is given in Section V. 

II. SVD and Matrix Sample Paths 

In this section, the SVD will be reformulated to allow 
greater freedom when computing the factors of a random 
matrix. Next, it will be shown how the SVD of a MIMO 
channel matrix is used to form a spatial filter that is matched 
to the radio propagation paths between transmitter and receiver 
arrays. The problem of tangled sample paths will then be 
introduced. 

A. The SVD Revisited 

Consider the set { £ C MxN \ k = 1, . . . , K). Each 
admits a SVD according to, 

ff( k ) = f/( fe )5]( fe )y( fe ) t (i) 

where U^ k) e C AIxM and V {k) e C NxN are unitary matrices 
whose columns, , ■ ■ ■ , u'y 1 j and , ■ ■ ■ , vffi |, are 

the left and right singular vectors of H^ k \ The matrix, 
£( fe ) 6 R MxN , is populated by the singular values of 



H (k) 



, H fe) , 



min(M,JV) 



along the diagonal and zeros 



elsewhere. Two strict identification conditions are imposed on 
the SVD: 

i. The singular values are real and non-negative. 

ii. The singular values are ordered along the diagonal as, 



■■> a 



(fe) 

min(M,A0 



(2) 



and the singular vectors are ordered accordingly. 

These strict identification conditions may be relaxed as fol- 
lows: 

(a) It is permissible to multiply by cxp (jOjji^ and 
t^ fe ) by cxp (jOyl^j without violating so long as 

4 k) is multiplied by exp (-j + 0$)). 

(b) The singular values may be permuted on the diagonal 
of £w without violating Q so long as the columns of 
t/W and are permuted accordingly. 

In what follows, the case of square matrices (M = N) may 
be considered without loss of generality (the generalization 
to rectangular matrices is easily accommodated). Points (a) 
and (b) can be incorporated in the SVD definition. Let 
p( k ) g ^NxN k e a unitajy permutation matrix such that post 
multiplication of a matrix by p( fc ' causes the reordering of 
the columns of the matrix. Let 6 fcNxN ^ e a diagonal 
unitary matrix given by, 



9 (fc) = diag 



e 3 N 



(3) 



where 9^ has its domain on [0, 2ir). By reforming the SVD 
factors as, 



the SVD becomes, 



U(k) _ rr(fe) y(k) v (k) H 

n — u ep Zj ep y ep ■ 



(4) 
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There are therefore an infinite number of valid decompositions 
for each 

B. SVD Beamforming 

The noise free MIMO system is modeled as, 



= (V (fc) E (fe V (fc)ff ) x (k) 



(6) 



where (x,y) € C JVxl are the transmitted and received signal 
vectors respectively. A treatment of the noisy case is beyond 
the scope of this paper and is not necessary for conveying the 
ideas presented hereafter. A MIMO channel can be resolved 
into a set of parallel independent channels if the transmitted 




Fig. 1 . Array beam-pairs associated with the two dominant singular channels 
of a 10 X 10 MIMO channel. Beam gain-patterns are plotted on a linear scale. 



symbols are spatially multiplexed and the received symbols 
are spatially demultiplexed as follows [1], 

-(ft) = jj(k) H H (k) (y(k)g(k)^j 

= EW£ (fc) . (7) 

Furthermore, spectral efficiency is maximized by waterfilling 
transmit power across the elements of a?w according the 
channel's singular values. Such a scheme is accomplished 
by computing the SVD of at the receiver and feeding 
back E( fe ) and to the transmitter before performing the 
waterfilled mux/demux communication of x^ k \ The feedback 
rate must be fast enough so that this can be accomplished 
before the channel has changed significantly. On the other 
hand, feedback requires use of the channel (overhead) which 
results in a loss of spectral efficiency. The optimal feedback 
rate is therefore a function of the singular values' temporal 
coherence and feedback overhead. A similar situation arises in 
beamforming for MIMO OFDM systems where the singular 
values of each subcarrier channel vary across frequency and 
optimal pilot density is a function of coherence bandwidth and 
pilot overhead. The mux/demux operations have an important 
beamforming interpretation [6]. Consider the simulated MIMO 
propagation scenario shown in Fig. Q. Two circular arrays 
with ten elements and 1A radii are set up in a simple Line 
of Sight (LoS) scenario with four scatterers. Propagation 
paths are discovered by raytracing between each pair of 
transmit and receive antenna. A simple path loss model is 
then applied to each path and the MIMO channel matrix 
and its SVD are computed. Pairs of left and right singular 
vectors are then used as weighting vectors to beamform at 
the transmitter and receiver [7]. The beam patterns associated 
with the two strongest singular values are shown. The first 
beam-pair focuses energy along the LoS path while the second 
beam-pair uses the scattering paths. The remaining beam- 
pair patterns (not shown) are unique and use the available 
multipath to provide other independent channels. Thus, the 
SVD relates the MIMO channel matrix to the underlying 
physical propagation environment by way of the spatial filters 
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formed by the singular vector pairs. It is important to recognize 
that a common phase rotation of the elements of a singular 
vector does not alter the associated beam's gain-pattern. Thus, 
even though there are infinite valid factorizations of a MIMO 
channel according to (0, there is a unique set of beam gain- 
patterns for a given propagation scenario. 

C. Tangled Sample Paths 

Consider a series, fft 1 -*) = {H^\ . . . ,H (A ' ) }, that is a 
sample path taken from a wide sense stationary and bandlim- 
ited random matrix process at a rate satisfying the Nyquist- 
Shannon sampling criteria. Such a sample path arises when 
sampling a MIMO radio channel across time or frequency. 
The SVD may be computed for each to yield a sample 
path decomposition, 

H {1...K) ^ u {l...K)^(X...K)^y{X...K) (g) 

where the sample paths of the SVD factors are, 

^•••*) = {[/«...,[/(*)}, 

E^ 1 -*) = {£«...,£(*)}, 

F (i...if) = | F (i) ) _ .. )V W}, (9) 

and sample paths of the singular channels are denoted by the 
triplet, 

,~ -*.(\...K) 



for i = 1, . . . ,min (M, N) 



(10) 



A SISO channel sample path, h^"' K , has a smooth evo- 
lution that is governed by the covariance of the underlying 
random process [8]. This random process arises from the 
physical process of multipath superposition. Since both SISO 
channels and MIMO singular channels arise from similar 
physical processes, it is expected that the singular channel 
sample paths should have a smooth evolution. However, when 
the strict identification conditions, (i) and (ii), are imposed 
independently from one sample to the next when computing 
(|8}, the resulting singular channel sample paths will be a 
tangled ordering of the natural singular channel sample paths. 
This will be made clear by example in subsequent sections. 
The tangled singular channel sample paths will be seen to not 
have a smooth evolution and will have an auto-covariance and 
a cross-covariance that are not consistent with the underlying 
random process. This tangling is a serious problem for closed 
loop MIMO communications because the feedback rate (pilot 
density) needed for channel tracking (channel estimation) is 
greater than what is inherently necessary. Consequently, the 
system's performance is diminished. 

III. Untangling Singular Channels 

A solution to the problem of tangled singular channel 
sample paths can be formulated by using the SVD with 
relaxed identification conditions to compute singular channel 
sample paths whose covariance is consistent with that of 
the underlying random process. The SVD formulation in l|5} 
requires a search over the space of permutation and rotation 



matrices for the factorization that is most consistent with those 
of the adjacent matrices' SVDs. Since it is the singular vectors 
that contain information about the propagation multipath, they 
will play the central role in finding the optimal factorization. 
As the search proceeds from one sample to the next, the 
singular channel sample paths will untangle from the muddled 
sample paths computed with the strict SVD identification 
conditions. 

A. Correlation Recovery Strategy 

A reference matrix, flW, and a target matrix, H^ k > , are cho- 
sen from the sample path, H^"" K \ The assignment = 
may be made without loss of generality. The SVDs 
of H<- R *> and i?W are then computed subject to the strict 
identification conditions. The triplet (f/ (R) , S (R) , F (R) ) will 
serve as the template to which the triplet (U^, Y,( k \ V^) is 
matched as closely as possible using permutations and phase 
rotations according to ©. Let the difference matrices between 
the singular vectors of the and H^> be, 



D 



a-) 



(fc) 
BP 



D 



(k) 



V^-V @ p (11) 



The matching can be formulated as a search over all possible 
triplets &y\ P' fc )^ with the objective of minimizing 

the sum of squared differences, 



mm tr 

p( fc ), e^',e' fc) 



(fcr n ( fe ) 



D\J' D 



u 



(k) H n (fe) 



D y y' D 



V 



(12) 



where the sum of squared differences can be expanded as, 
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Since terms one, three, four and six on the right side of Jl 3I > 
are not affected by minimizing the sum 

of squared differences is equivalent to maximizing the sum of 
singular vector correlations (terms two and five), 



max TZe 



R 



R 



(fc) 



where the correlation matrices, 

= U^ H U e p , R (k) = V^ H V eP , 



(14) 



(15) 



have elements that comprise the set of all possible inner 
products of the reference and target singular vectors. The 
maximization in (114-i implies a search over the infinite set 

' ci( fc ) p(ky 



{( 



non, [eP,eP,P(v 



Fortunately, the global optimum solu- 
, can be found by searching over a 
are known. As will 



finite set if either P< fe ) or (e^ , 9^ 
be shown, an elegant search strategy is to find P 1 -^ assuming 
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that an appropriate ,Qy^ exists and then to compute 

(e^.G^) analytically given pW. 
The search for p( k > proceeds as follows. A simple inter- 

(k) 

change columns i and j of Uq P is equivalent to interchanging 
columns i and j of R^f' 
Consequently, the measure of improvement in, 



(and similarly for V@p and Ry). 
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that results from an interchange of columns i and j is, 
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Phase rotations of the singular vectors according to @ by the 
guarantee that the greatest possible 



appropriate 
improvement, 
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is achievable, thus permitting the search for P^ to be 
executed by column swaps without jointly searching for 6 
and Qy . The maximization problem is then, 



max tr 



R 



R 



(fc) 



(19) 



The search for can be performed iteratively by column 
swaps using ( I18> as a metric for directing the search. The 
details of an efficient algorithm for implementing this search 
strategy are given in Section III-B. Once p( h > has been found, 
finding the associated (q^,Q^J is simply a matter of 

computing the angle between singular vectors of and 
the permuted singular vectors of H 1 -^ according to, 



Jk) 



u) u) 



1, 



(20) 



and likewise for {0y ]} where arg (•) = 1m {In (•)} yields the 
radian angle of a complex scalar. These left and right phase 
rotation matrices are then applied to T,^ as required by l|5} 
to compensate for the conjugate phase rotation being applied 
to the corresponding singular vectors. The resulting match of 
'U^l,V^) to ([/( R ) , yW) is optimal in a least squares 



sense. 



Thus far, the discussion has assumed noise-free sample 
paths. The case of an additive noise model, y w = H^x^ k ' + 
n^ k \ is relevant and can be addressed in part by using the 
weighted singular vectors, 



jj(k) = Jj(fc) ^(.k)y /2 ^ y-(fc) = y(k) /^(fc)^ 



1/2 



(21) 



Q(fc) p(k) \ rp n j s biases me search 
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in favor of those singular vectors with greatest signal to noise 



Table l : Untangling Algorithm 
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Seed the reference: H<- R > = H il] 
Compute: H m --5 y(R) s (E)y(R)* st a W > . 
for/: = 2,3,..., K 
Compute: #<*> 5-5 ErWEMvM', s.t of 3 > 



> a 



Set PV" = I 
Set improving = /nie 
while improving — true 
Set improving = false 
for i = {1,2, ...,JV} 
SetAfS (i.i) = 
tm j = {!,%,. ..,N}^i 
Compute A<*j(i, j) 
if A (k) (i i) > A (tl) 



jfjM = y(R)"y(t) 



Set A 



(*>j)l jmax = j 



M 

fj il = A" 

if A^ > 

Interchange columns i and j max of J% , JjL , P'*' 
Set improving = true 
Compute: {^Ji = 1, ...,n\ and {e^Ji = 1, ...,Jv] 
Form: 0*' and ef J 
Compute: E^, s£ 
Set !J( R ) = [/<*> E<*> = £ ( *'. y' R 



p 



ratio. This is not all that can be done to counteract the noise. 
However, a proper treatment of noisy sample paths is beyond 
the scope of this paper. 

B. Untangling Algorithm 

The algorithm given in Table 1 untangles the singular 
channel sample paths from the sample paths computed with 
the strict identification conditions. The algorithm's core (lines 
4-21) is an efficient implementation of the correlation recovery 
strategy discussed in Sections III-A. The number of possible 
permutations on the singular vector ordering is Nl. However, 
the core algorithm guarantees that the least squares solution 
for P^ is found in at most N (N — 1) column swaps. The 
essential aspects of the core algorithm are: 

i. Column i (pivot column) is swapped with the column 
that would result in the maximum improvement to i ll SL 
whenever an improvement is possible. 

ii. Each column becomes the pivot column in cyclic order 
and the algorithm proceeds iteratively until no further 
improvement can be made. 

More generally, the core algorithm can be applied to any 
problem where the trace of a matrix must be maximized 
under column (or row) permutations. The algorithm's shell 
(lines 1-3,22) parses the sample path, H^"" K \ using a sliding 
reference, H^> = H^ k ^ x \ The SVD factors computed by the 
core algorithm for each matrix sample are therefore dependent 
on the reference seed-matrix's factorization (line 2). This seed 
matrix has its SVD computed according to strict identification 
conditions on ordering and phase. Any ordering and phase for 
the seed reference's SVD would also suffice. Whatever the 
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initial order and phase, they determine the relative ordering 
and phase-offset of the untangled singular channel sample 
paths. As will be discussed in Section IV, the untangling 
solution is unique when conditioned on the seed matrix's 
singular value ordering and phases. 

The Untangling Algorithm's performance can be poor if 
either of the following two situations arise and are not treated 
appropriately. The first situation is that of a matrix with close 
singular values. In this case, the SVD algorithm becomes 
unstable and slight perturbations of the matrix cause wildly 
different decompositions. The second situation (an extreme 
case of the first) is that of a matrix with repeated singular 
values. In this case, the singular vectors associated with the 
repeated singular values are indeterminate and contain no 
information. The case of equal singular values is very unlikely 
in practical applications where the H^ k > are sampled from 
a random processes occurring in the physical world. These 
two situations may be handled by implementing the following 
simple modifications to the Untangling Algorithm: 

i. Set a threshold, ^S k \ for the tolerable difference between 
singular values returned by the SVD algorithm. For 
example, 

min (a[ k) - cr^) < e, i = 1, . . . , N - 1 (22) 
where e is a small positive number and a[ k ^ > • • • > 

N ■ 

ii. If the threshold is violated then attempt a match but 
display a warning message and blacklist the offending 
matrix so that it cannot be used as a reference. 

It is possible, though unlikely, that a tie is computed for 
different column swaps' improvement metrics. This situation 
is averted in line 14 by the greater-than logic which causes 
the first of the column swaps ties to be used. Even so, a tie 
should be treated in a manner similar to the case of repeated 
singular values by giving a warning, attempting a match and 
blacklisting the offending matrix. 

IV. Results 

In this section, an example is given to demonstrate the effec- 
tiveness of the Untangling Algorithm in computing smoothly 
evolving singular channels from a fading MIMO radio channel. 
A comparison is made between tangled and untangled sample 
paths. The covariance of the untangled singular value sample 
path is compared with those of the SISO channel sample paths 
and tangled singular value sample paths. Joint and marginal 
distributions of the untangled singular values are compared 
with those known from random matrix theory. 

A. MIMO Channel Synthesis 

Consider the case of a MIMO radio distortion process 
between a fixed transmitter and a mobile receiver. It may 
be assumed that the random MIMO distortion process, H, 
is a time varying matrix of i.i.d. zero mean Gaussian random 



variables, fry. The temporal auto-covariance of each SISO 
channel realized from H may be given by [8], 

c(h%,h^ T) )=Jo(2nf d T) (23) 

where Jo ( x ) is the Bessel function of the first kind and 
Id = v/\ c is the maximum Doppler frequency resulting from 
a receiver moving with velocity, v, and a carrier frequency 
with wavelength, A c . The assumption that the SISO channels 
are independent and identically distributed is justifiable in 
many communication scenarios. The temporal auto-covariance 
given in d23i arises from non-line-of-sight scenarios with 
rich multipath propagation. This bandlimited random matrix 
process may be sampled over time to yield a sample path, 
Jj0~-K) Realizations of h^" can readily be generated in 
computer simulation by a multipath sum-of-sinusoids model 
[9]. This model does not generate MIMO channels directly 
from antenna array geometries and a scattering topology as in 
the case of Fig.^ However, array geometries and a scattering 
topology are associated with the model's output by way of the 
the random variables used to seed the paths' phase and arrival 
angle. For the analysis that follows, MIMO channel sample 
paths were generated according to the sum-of-sinusoids model 
with these parameters: 500 paths, fd = 15 Hz, A c = 0.125 m, 
and T s = 3.3 ms. 

B. Untangled Sample Paths 

Fig. |3 a) shows the singular sample paths for a realization 
of H^ 1 " ' where M = N = 3. SVD's were computed 
according to strict identification conditions on ordering and 
phase with the result that singular value sample paths form 
separate layers. Fig. |2jb) shows the average singular vector 
correlation between adjacent samples, 

R (k) = ±tr (ir« H E#] + VW'VM) . (24) 

The correlation profile indicates that although the singular 
vectors of the k th and (k — l) th samples are initially corre- 
lated, the strict identification conditions soon cause a loss 
of correlation. Fig. |3a) shows the magnitude of the singular 
values after applying the Untangling Algorithm to the same 
MIMO sample path as was used for Fig. |2ja). Although the 
singular values magnitudes have not changed, the traces that 
connect them are very different than what they were before 
untangling. The singular value sample paths are now observed 
to weave. When the traces cross, they do so in natural smooth 
transitions. The correlation profile in Fig. |3jb) shows that the 
untangled singular vectors maintain high correlation with those 
of the earlier sample. An important aspect of the untangling 
process is that the singular values are allowed to become 
complex valued. Fig. |4] shows the untangled singular value 
sample paths as a function of time. These traces reflect a 
natural smooth evolution in accordance with the bandlimited 
underlying random process. 

Furthermore, the untangling solution is unique when condi- 
tioned on the seed matrix's singular value ordering and phases. 
For example, consider the forward untangled singular value 
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Fig. 2. (a) Singular value sample paths of a spatially white 3x3 MIMO 
channel computed using strict identification conditions on singular value 
ordering and phase, (b) Average singular vector correlation between adjacent 
samples. 
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Fig. 3. (a) Magnitude of the untangled singular value sample paths of a 
spatially white 3x3 MIMO channel, (b) Average singular vector correlation 
between adjacent samples. 




Fig. 4. Smooth evolution of the untangled singular value sample paths from 
a spatially white 3x3 MIMO channel. 



sample path £(^> 2 >— ^0 generated by parsing H^ 1,2, '" K ^ and 
the reverse untangled singular value sample path, £C^ 2 >— K > 
generated by parsing H^ K ' ' 2,1 \ It was verified that, 

E a.a,...Jf) = Pr q ^(5,. .,2,i) q ^ Pr (25) 

without exception (except for numerical precision errors). Here 
r is a diagonal unitary matrix that derotates the singular 
values in Tr^> so they are real and non-negative, P r is the 
permutation matrix that is then used to perform the reverse 
one-to-one mapping of the singular values in £^P to those 
in S£!y and represents element-wise matrix multiplication 
over the entire sample path. The existence of a unique solution 
for a given MIMO sample path is consistent with the specific 
underlying physical process giving rise to that sample path. 

C. Sample Path Covariance 

Fig- compares the magnitude of the temporal covariance 
of a SISO channel sample path, a raw singular value sample 
path, and an untangled singular value sample path. These are 
computed from 3x3 sample paths of length K = 10 000. In 
general, each of the min (M, N) raw singular value sample 
paths has a slightly different covariance function whereas the 
covariance of each of the untangled singular value sample path 
is the same. For simplicity, the covariance shown for the raw 
singular values is the mean of all covariances. 

Figure |2a) shows that the SISO channel auto-covariance 
matches that expected from d23l > and reaches 0.7 at t = 
0.012 s. The auto-covariance for the raw singular value sample 
path decays slightly faster than that of the raw singular 
channels and reaches 0.7 at r = 0.006 s. It is remarkable 
to observe that the untangled singular channel's covariance 
decays far less rapidly than that of the SISO channels and 
reaches 0.7 at r = 0.035 s. The untangling process has 'slowed 
down' the perceived channel dramatically. This result implies 
that channel tracking (estimation) can be done with a lower 
feedback rate (pilot density) when using the untangled singular 
channels instead of the raw MIMO channels directly. This 
promise of improved efficiency can be realized by implement- 
ing the Untangling Algorithm at the receiver. Furthermore, 
a transmitter with a priori knowledge of the singular chan- 
nels' covariance matrix function has an extended horizon for 
predicting the channel. This is a benefit because latency due 
to channel estimation and feedback carries the consequence 
that the transmit beamforming coefficients may no longer be 
matched to the current channel. The transmitter may then use 
the most recent untangled channel state information to predict 
the current channel state for beamforming and waterfilling. 

Figure |3b) compares the magnitude of the temporal cross- 
covariance of a SISO channel sample path, a raw singular 
value sample path, and an untangled singular value sample 
path. As expected, the SISO channel sample paths are uncor- 
rected. The tangled nature of the raw singular value sample 
paths manifests itself as a non-zero covariance at lag r = 0. 
The fact that the untangled singular value sample paths have 
no covariance verifies that the untangling process has resolved 
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Fig. 5. Covariance functions computed from a spatially white 3x3 MIMO 
channel (a) auto-covariance (b) cross-covariance. 



sample paths whose evolutions are independent of each other. 

D. Sample Path Density Functions 

Analytic results are known for the distributions of ordered 
and unordered eigenvalues of the Wishart matrices [1], [10]. 
This is relevant to the topic of MIMO singular channels 
because W^ k > = H^H^ is drawn from the Wishart 
distributed random process W = HH H . The eigenvalues 



of W^ k \ 



. A 



M 



are simply the magnitude squared 

of the singular values of H<- k \ Given that e C MxN , the 
joint density function for the unordered eigenvalues of W is, 



M 



\N-M 

,A i .) = M exp(-E^)n (M _; ) , (iV _. )! 
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(26) 

while the marginal density function of the unordered eigen- 
values of W is, 
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Tr E 77" 
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(27) 



where L k (x) is the associated Laguerre polynomial with 
Rodrigues representation, 



LI (x 



) = E(-!) 



(n + k)\ 



(n — m)l (k + m)! m! 



(28) 



Fig. [6] shows the joint density function computed from 
(I26> for N = M = 2. Fig. shows the joint distribution 
of the untangled eigenvalues for a realization of ffi 1 - 10000 ) 
and computed using a histogram over the same domain as 
the function in Fig. [6] The untangled distribution is matches 
the theoretic density function and is completely different 
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Fig. 6. Joint density function given by 1261 for the unordered eigenvalues 
of W where H G C 2x2 and H ~ M (0, 1) 




Fig. 7. Joint distribution of the eigenvalues of from an untangled 

sample path H^ 1 - 1 ^ where ff< fe ) e c 2x2 and ~ AT (0,1). 



from the joint distribution of the ordered eigenvalues (not 
shown). It is difficult to show the joint density function for 
min (M, N) > 2 but the marginal density function for these 
higher order cases can be shown as is done in Fig. [8] for 
M = 2 and N = 3. Again, there is very good agreement 
between the untangled eigenvalue distribution and the theoretic 
distribution. It is expected that the untangled eigenvalues 
should have distributions agreeing with those given in l !26> and 
l !27t because the untangling operation recovers the ordering 
randomness of the underlying i.i.d. Gaussian process. 

V. Conclusion 

This work has addressed the problem of how to compute 
the SVD's of a random matrix sample path in a way that 
preserves the covariance of the underlying random process. 
The solution to this problem was developed as follows: (i) 
a SVD with relaxed identification conditions was proposed, 
(ii) an approach was formulated for computing the SVD's 
of two adjacent matrices in the sample path that maximizes 
the correlation between corresponding singular vectors of the 
two matrices, and (iii) an efficient algorithm was given for 
untangling the singular value sample paths. 



default 




[9] R. H. Clarke, "A statistical theory of mobile-radio reception," The Bell 
System Technical Journal, vol. 47, pp. 957 - 1000, July - August 1968. 
[10] A. M. Tulino and S. Verdu, Random Matrix Theory and Wireless 
Communications. Now Publishers, 2004. 



Fig. 8. (i) Marginal density function given by 1271 for unordered eigenvalues 
of W where ff £ C 2x3 and H ~ A^(0, /) and (ii) computed eigenvalue 
marginal distribution for W^°' from an untangled sample path H0---K) 
where £ C 2x3 and ifW ~ Af(0,I). 



The algorithm's effectiveness was demonstrated on i.i.d. 
Gaussian MIMO channels. It was shown that the algorithm 
resolves smoothly evolving singular channel sample paths that 
are in accord with the stochastic structure of the underlying 
random process. Furthermore, the algorithm gives a unique so- 
lution conditioned on the seed matrix's singular value ordering 
and phases. 

It was shown that the untangling process dramatically 
increases the coherence period (or bandwidth) of the singular 
channels. A primary application of the algorithm is in MIMO 
radio systems. The benefit promised by using SVD untangling 
in these systems is that the fading rate of the channel's SVD 
factors is greatly reduced so that the performance of channel 
estimation, channel feedback and channel prediction can be 
increased. 
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